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The dynamical properties of a recently introduced phenomenological model for high temperature 
superconductors are investigated. In the clean limit, it was observed that none of the homogeneous 
or striped states that are induced by the model at low temperatures can reproduce the recent angle- 
resolved photoemission results for LSCO (Yoshida et ai, Phys. Rev. Lett., 91, 027001 (2003)), 
that show a signal with two branches in the underdoped regime. On the other hand, upon including 
quenched disorder in the model and breaking the homogeneous state into "patches" that are locally 
either superconducting or antiferromagnetic, the two-branch spectra can be reproduced. In this 
picture, the nodal regions are caused by d-wave superconducting clusters. Studying the density 
of states (DOS), a pseudogap is observed, caused by the mixture of the gapped antiferromagnetic 
state and a d-wave superconductor. The local DOS can be interpreted using a mixed-phase picture, 
similarly as observed in tunneling experiments. It is concluded that a simple phenomenological 
model for cuprates can capture many of the features observed in the underdoped regime of these 
materials. 
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I. INTRODUCTION 



Since the discovery of the high-T c cuprates it has been 
widely suspected that key information for the under- 
standing of superconductivity lies in the curious regime 
located between the insulating (antiferromagnetic, AF) 
and metallic (superconducting, SC) phases. Whereas 
these two phases are individually fairly well understood 
phenomenologically, the intermediate regime is, even af- 
ter several years of research, still to a large degree mys- 
terious. Experiments have often only been able to de- 
cide what this state is not - for example, it is not a sim- 
ple Fermi liquid or any other clearly defined, well-known 
state. This has given rise to a variety of proposals in 
describing this peculiar regime, often in terms of "ex- 
otic" order scenarios such as competing charge-density 
wavei (CDW) or pair-density wave states^, staggered- 
fiux phases 3 -, spin-Peicrls4 states and, among others, or- 
bital currents^. The much-debated pseudogap (PG) in 
the DOS that appears in this phase may then be regarded 
as a manifestation of a hidden order. Alternatively, it has 
been suggested that the strange phase interpolating be- 
tween AF and SC states may be characterized by a partic- 
ularly strong attraction between charge carriers, strong 
enough to drive T c down, and leading to a state of pre- 
formed, but yet uncondensed, pairsi. In this scenario, the 
PG temperature Tpq signals the onset of pairing fluctu- 
ations. On the other hand, from the experimental view- 
point the regime between the AF and SC phases is often 
described as "glassy", with slow dynamics^. The long 
discussions on these issues show that an understanding 
of the low hole-doped cuprates has not been reached yet, 



In recent years, progress in both sample preparation 
and measurement techniques has led to very interesting 
experimental results for this fascinating phase. Espe- 
cially important in this context are the insights gained via 
angular resolved photoemission spectroscopy (ARPES) 9 , 
since this technique allows for a direct tracking of the 
Fermi surface (FS) - if it exists - and therefore provides 
crucial information for the evolution of the metallic phase 
from its insulating parent compound. Moreover, in the 
particular case of Lai_ x Sr x Cu04 (LSCO), a relatively 
simple single-layer cuprate, those data have been ob- 
tained over the whole doping range x, starting from the 
antiferromagnetic Mott insulator at x=Q up to the "op- 
timal" doping x=0.15ia. These new results apparently 
show that portions of a FS (nodal regions) are present 
even for the smallest doping levels considered, such as 
x=0.03, where the material is not superconducting™. 
This suggests that even the slightly doped Mott insulator 
is emphatically different from its half-filled parent com- 
pound, and that it should be best understood as some 
sort of insulator with coexisting (embedded) metallic re- 
gions. In fact, the doping evolution of ARPES data- 
displays two different branches, one evolving from the 
insulator deep in energy and the second created by hole 
doping and containing nodal quasiparticles. In addition, 
it was found that the FS fragments, which are located 
around (ir/2,n/2), do not expand as the hole density is 
increased - as one would expect for a conventional metal- 
lic state - but rather acquire more spectral weight. Only if 
the temperature is increased does this FS arc widen, rem- 
iniscent to the closing of the gap in superconductors. As 
it will be argued in this paper, the most natural explana- 
tion for all these results is that there arc already d-wave 
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paired quasiparticlcs present in the glassy phase. The 
simplest such picture is one of phase-separation (PS), 
where the metallic and insulating states inhabit spatially 
separated (nanoscalc) regions, a view recently introduced 
by Alvarez et alH. This explanation is against the ex- 
otic homogeneous states proposed as precursors of the 
SC state, but nevertheless the reader should note that a 
mixture AF+SC is highly nontrivial as well, and in many 
respects it is as exotic as those proposals. For instance, 
this state has "giant" effects, as recently discussed^, re- 
sults compatible with those found experimentally by Bo- 
zovic et alm& and othersi^i^ in the context of the giant 
proximity effect. 

One might also note that similar concepts have been 
proposed in the context of the metal- insulator (M-I) tran- 
sition in the manganites, and have found widespread 
acceptance ther o 16 i 17 . For this reason, there exists 
a concrete possibility to frame the M-l transition in 
strongly correlated systems in a unifying picture of phase- 
separation and subsequent percolation. Indeed, such a 
proposal has been madei& (by some of the authors), 
claiming that such a mixed-phase state is the conse- 
quence of impurity effects acting in regions of the phase 
diagram, where the AF and SC phase are very close 
in energy. In fact, quenched disorder has such an ex- 
traordinary and unusual influence because of the pos- 
sible first-order nature of the clean limit M-I transition 
(which is not experimentally accessible), and the effective 
low-dimensionality of the transition-metal oxides under 
consideration. Note that in some scenarios^^*^ a similar 
sharp transition AF-SC is invoked, but the influence of 
quenched disorder has not been explored. In our case, 
quenched disorder is crucial to produce the intermediate 
state between the AF and SC phases. 

An important experimental tool in detecting coex- 
isting phases is scanning tunneling microscopy (STM), 
which provided clear indication for mixed-phase states 
in manganites years agoi&ii. More recently, similar 
results have also been available in slightly underdoped 
Bi 2 Sr2CaCu208+,5 (BSCCO)2i, suggesting the existence 
of metallic and insulating phases on nanoscale (typi- 
cally 2-3 nm) regions. Unfortunately, it is not possible 
to extend these measurements into the strongly under- 
doped regime, as BSCCO becomes chemically unstable. 
However, such STM results were recently reported by 
Kohsaka et al~ for Ca2_ x Na :E Cu02Cl2 in the neigh- 
borhood of the superconductor-insulator transition point 
at :r=0.08. The compound Ca 2 - x Na :i: Cu02Ci2 is struc- 
turally closely related to LSCO, and therefore offers a 
unique opportunity to compare real- and momentum- 
space imaging. Kohsaka et al. identified separate metal- 
lic (SC) and insulating areas of approximately 2nm in 
diameter and demonstrated that the two phases in ques- 
tion are only distinguished by a different ratio of metallic 
and insulating clusters. Together with the BSCCO data 
this points toward a mixed-phase description of high tem- 
perature superconductors. 
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challenged since they are surface probes and not bulk in- 
vestigations. However, recent bulk measurements using 
Raman scattering have provided results that also may 
best be interpreted in terms of mixed-phase states and 
have given further support for such scenarios. Machtoub 
et alm& have studied LSCO at x=0.12 in the SC phase and 
in the presence of a magnetic field. The results were inter- 
preted in terms of an electronically inhomogeneous state 
in which the field enhances the volume fraction of a phase 
with local AF order at the expense of the SC phase. The 
same conclusions were reached upon studying neutron 
sacttering in underdoped (a;=0.10) LSCO, where it was 
claimed that the observed appearance of an AF phase as 
the SC is suppressed by an applied magnetic field points 
towards coexistence of those two phases^. Recent in- 
frared experiments on the Joscphson plasma resonance in 
LSCO have also suggested a spatially inhomogeneous SC 
stated. The same conclusion was reached investigating a 
high field magnetoresistance2&. Studies by Keren et aim! 
using fj,SR techniques have also led to microscopic phase 
separation, involving hole-rich and hole-poor regions. 

Clearly, the notion of electronic phase separation^ has 
a long history in the cuprate literature, going back at 
least to the original proposal of the stripe state^ffiiSiiS. 
Although stripes as originally envisioned should only ap- 
pear for very specific doping levels, this may simply be a 
consequence of the approximation used and one could 
expect stripe-like correlations or corresponding charge 
fluctuations to appear for a broad array of doping frac- 
tions. Although stripes are only one possible avenue to 
introduce mixed-phase tendencies, they have nevertheless 
attracted considerable interest. This is understandable 
since in the mid 90's experimental evidence from inelas- 
tic neutron scattering seemed to demonstrate stripe-like 
charge-order in Nd-doped LSCO at hole doping x=l /8~. 
However, the more recent investigations mentioned here 
suggest a broader picture of phase separation, with nan- 
ocluster "patches" of random sizes and scales, rather than 
more organized quasi-one-dimcnsional stripes. 

Even if a (general) mixed-phase scenario is accepted 
owing to experimental evidence, it is still not clear 
whether or not the electronic inhomogencity is caused 
by the (supposed) inherent tendency of strongly-coupled 
fermionic systems to phase-separate or by the imperfect 
screening of the dopant ions due to the proximity of the 
Mott insulating phase and concomitant strong disorder 
potentials. Those two scenarios are not mutually exclu- 
sive, but put different emphasis on the aspect of random 
(chemical) disorder: whereas in the latter scenario it is 
thought that disorder is strong enough to overcome the 
tendency of fermionic ensembles to form a homogeneous 
system, particularly when phases compete, in the first 
one chemical disorder may act as a mere catalyst of al- 
ready present tendencies, possibly leading to a pinning of 
stripes/charge-depleted regions. If, as one might assume, 
underdoped cuprates lie in between those extremas, some 
materials could be more influenced by disorder than oth- 
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systems when studying lightly-doped transition metal ox- 
ides. 

Our goal here is to investigate the spectroscopic prop- 
erties of systems with competing AF and SC order, and to 
compare them to recent ARPES and STM experiments, 
in order to develop a coherent understanding of cuprates 
from the mixed-phase scenario point of view. The present 
paper builds upon the recent publication by Alvarez et 
alH, where the proposed state was described in detail 
and "colossal" effects for cuprates were predicted. The 
approach we followed assumes as an experimental fact 
that AF and SC phases compete, and addresses how the 
interpolation from one phase to the other occurs with in- 
creasing hole doping, once sources of disorder are intro- 
duced. It is concluded here that recent LSCO ARPES 
results can be neatly explained by the mixed-phase sce- 
nario. The paper is organized as follows: In Section II 
the formal aspects of the problem are described, includ- 
ing models and approximations. Section III addresses the 
ARPES data for (i) a variety of uniform states (showing 
that none fits the experimental results) and (ii) in the 
presence of quenched disorder that leads to the mixed- 
phase state. Results for the latter are found to reproduce 
experiments fairly well. In Section IV, our results for the 
local DOS are presented and discussed. Conclusions are 
provided in Section V. 



II. THEORY AND METHODS 

A minimum model for describing the interplay between 
AF, SC (and possibly charge-order (CO)) can be found 
in the extended Hubbard model -Htjv, defined as 



H 



UY 



<1J>,CT 



(l) 



where c- la are fermionic operators on a two-dimensional 
(2D) lattice with N=LxL sites, t is the hopping between 
nearest-neighbor (n.n.) sites i, j, and serves as the energy 
unit. Ui is the usual Hubbard term and Vy>0 describes 
an effective n.n. attraction. The local particle density 

^i = 5Z<T C L c i(T ' s regulated by the chemical potential m, 
which, like the interactions Ui and Vy, is allowed to vary 
spatially. To address the properties of this model, the 
four-fermion terms in i?uv are subjected to a Hartree- 
Fock decomposition, 

n iT n u ->■ ra iT (ra u ) + (ra iT )ra u - (n^)(mi), 

m«j A ij c ji4r + A ij c iT c ji " l A ijl 2 > ( 2 ) 

where we have assumed that the predominant tendency 
in the Vy term is towards superconductivity rather than 
particle-hole pairing, which would favor charge-ordering 



new (site dependent) order parameters (o.p.), m\ and 
A„: 



mi = (n iT ) - (mi) 
A ij = (ciTCji)- 



(3) 



nii is the local magnetization, and Ay is the SC 
o.p., defined on the link ij. The Hamiltonian 
i/uv=^HF=^HF+^HF now * s quadratic in electron op- 
erators and written as 



Hhf = -t ( c L c j CT +H.c.) - Y^W™ 

<ij>,er iff 



<ij> i 
+ 5>y|Ay| 2 + 1/4^(71?-™?), (4) 
<ij> i 

after we have introduced the local spin operator sf=n^- 
mi ■ Hhf is effectively a single-particle Hamiltonian with 
2N basis states, which can be readily diagonalized using 
library subroutines in the general case. It needs to be 
stressed that the third line in Eq.(0} (=H^ F ) contains 
c- numbers only, and no operator terms, unlike the first 
two rows (=i/{j F ). It is, however, the interplay between 
the second line in Eq.|0J, which tends to increase the 
o.p. amplitudes and the third one, which enforces an 
energy penalization for large values of |Ay|, m,\ that de- 
termines their actual values. We have resorted to two dif- 
ferent approaches in studying Hhf , a conventional mean- 
field method and a Monte Carlo (MC) technique, both 
of which we will describe in detail below. 



A. Mean-Field Method 

The mean-field method relies on the self-consistency 
condition Eq.@ to determine the appropriate values of 
m\ . Ai, and also n- la . Diagonalization itself amounts 
to performing a slightly modified Bogoliubov-de Gcnncs 
(BdG) transformatiorji^iSiiSi, where the electron opera- 
tors Ci CT are expressed in terms of new quasiparticles 7 n(T , 
defined as: 



n=N 



71=1 

C U = E'f 6 "^)^™^ +<+w( i )7 r l T }- (5) 



a„(i) and b n (i) in JSJ) are complex numbers and are cho- 
sen so that a Hamiltonian that is diagonal in 7 ra(T emerges. 
In this approach, the o.p. amplitudes are determined via 
the self-consistency condition © , which allows to express 
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Mi): 

N 

Ai J = 5Z a «( i ) a «+Jv(j)(7nT7 I t lT > - ^+jv(i)&n(j)(7ij.7n|), 

TV 

n n = 51 l a n(i)r(7l T 7„ T ) + I & n+A r ( i )| 2 (7„i7 1 t u ), 

n=l 
iV 

n U = Yl l fo «( i )| 2 (7li7„ T ) + |a„ + j V (i)| 2 (7„i7 r t li )- (6) 

n=l 

where (7i -7„ (J )={l + e' 35 '™''} -1 , i.e. the Fermi function 
/ with properties f(x)=l-f(-x), p=l/T the inverse tem- 
perature and E na are the eigenvalues of 10} . For the 
time being, however, we will work at T=0, which sim- 
plifies those relations considerably, and also guaranties 
better convergence of the self-consistent loop. 

The self-consistent set of Eqs. J3J), @, ©, and © can 
be solved in an iterative procedure, provided a meaning- 
ful starting wavefunction {a ra (i)}o is chosen. From the 
knowledge of the resulting cigcnfunctions, all relevant ob- 
servables such as the single particle spectral function 

A(k,u) = --ImG(k,L>), (7) 

where G(k, u) is the (retarded) Green's function and the 
local DOS 

iV(i,cj) = ImG(i,w) (8) 

can be calculated in a straightforward fashion. Those 
two observables are the crucial quantities with regards 
to ARPES and STM experiments, respectively. Although 
this approach is easy to implement and well-established, 
we want to mention here that it has its pitfalls, which 
lie in the correct choice of {a„(i)} . This is particu- 
larly problematic in the regime of small hole doping <5=1- 
"KfXa n ij <5<^1, which is known to lead to "stripe" states 
for not too small values of U, provided a sufficiently cor- 
rect selection of the seed functions is made. The exact 
nature of those stripe states, however, might depend on 
the linear lattice dimension L; also a possible interme- 
diate state between the undoped AF insulator and the 
stripe state is beyond the grasp of BdG (in the clean 
limit, at least). Of course, it is this particular regime that 
is of interest in this work - and for cuprates in general 
- and therefore we have in addition resorted to a novel, 
alternative approach, which does not suffer from the dis- 
advantages of the self-consistent technique, but provides 
an accurate, unbiased solution at any temperature. This 
is the MC technique described in the next subsection. 

B. Monte Carlo Procedure 

For this purpose, Hamiltonian (@J) is studied using a 
conceptually different approach, which stresses the im- 
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which only play a minor role in the self-consistent approx- 
imation. The MC technique also offers the additional ad- 
vantage of treating Ay as a complex variable, and allows 
us to write Ay = Ay |e 4 ^ ij , where </>y is the phase associ- 
ated with the bond between sites i, j. We will change 
our notation from here on and write Ay=Aio,, where a 
is a unit vector along the x or y directions. In a simi- 
lar fashion, we write <f>f instead of <^>y and also assume 
Vi a =Vi. Therefore, wc regard Ai Q as a site variable and, 
thus, finally write Ai Q . = |Ai|e 1< ^ i . 

Within the MC method, the partition function Zjjf 
pertaining to HhFj given as 



N r 1 „ r°° 

^hf = I / eW 3 ^™' / dlAile-^^IA.I 2 x 

i=i J- 1 •'0 

r dtfdflZcdm^UAil},^}), (9) 
Jo 



is calculated via a canonical MC integration over the am- 
plitudes mi, |Ai| and the phases In Eq.©, the 
C/ini-term from Eq.Q) has been absorbed into the lo- 
cal chemical potential. The purely electronic partition 
function Z c =Tr{e _,3ffHF } is obtained after diagonaliz- 
ing i?H F for a given, fixed set of those parameters and 
this diagonalization is responsible for limiting the lat- 
tice sizes. In fact, we chose another simplification and 
perform the "magnetic" integration over the signs only, 
rather than the amplitudes as well. This amounts to 
replacing the o.p. m\ by an Ising spin S?, which can 
take only the values ±1. We will also replace L^mi by a 
single term Jj, without loss of generality (i.e. a certain 
value of Ji corresponds to a certain value of U\ and vice 
versa). Compared to the standard procedure described 
above, the MC technique has some distinct advantages, 
in describing both the AF as well as the SC degrees of 
freedom (d.o.f.), since its results do not depend on the 
initial configuration. On the other hand, the MC tech- 
nique requires a very large number of diagonalizations, 
typically 10 6 rather than the ~10 2 iterations necessary 
to achieve self-consistency in BdG methods. This renders 
calculations on system sizes common in BdG (A^IOOO) 
impossible^i. 

The quantities of interest, such as the one-particle 
spectral function and the local DOS, can be evaluated 
in a straightforward fashion, cither directly or via the 
Green's function, which can be easily derived from the 
MC process. Similar techniques have been pursued in 
the double-exchange model, which is relevant for the 
manganites, and an in-depth description can be found 
in Refll6l. Further (technical) information with respect 
to the calculations presented here is provided in the ap- 
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FIG. 1: MC phase diagram for Eq.Q without disorder at low 
temperatures, using V=l-J/2, reproduced from Ref.ll2t Five 
regions were observed: AF, d-SC, stripes, coexisting SC+AF, 
coexisting stripes+SC, and metallic. Yellow dots indicate 
where A(k,a>) was calculated in the present work (see Fig|5J. 
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FIG. 2: A(k,io), evaluated via MC, on an 8x8 lattice 
for (a) (J,V,/i) = (0,1,-1) (SC state), (b) ( J, V, y) = 
(0.6,0.7,-0.4), coexisting AF/SC state, (c) (J, V, fi) = 
(0.7,0.65,-0.3) (AF), and (d) (J,V,fi) = (1,0.5,-1.2), 
striped state. 



III. SPECTRAL FUNCTIONS IN THE 
PRESENCE OF COMPETING STATES 

In this section, we will present the analysis of the one- 
particle spectral function for several regimes of the phase 
diagram of i?HF- For general doping and interaction val- 
ues this can only be done with the MC routine described 
above. 



A. Clean System 

The phase diagram of Hup for the clean case (i.e. 
without quenched disorder) was presented in Ref . and 
is reproduced for the benefit of the reader in Fig.Q] The 
figure shows two "paths" , which describe the transition 
from the AF to the SC phase. The first one crosses a 
region of long-range order with local AF/SC coexistence, 
whereas the second one involves an intermediate "stripe" 
state^i. We do not discuss here the exact nature of the 
stripe state, which may be horizontal or diagonal, de- 
pending on parameters such as doping and lattice size. 
For our purposes it is sufficient that an inhomogeneous 
state - stripe, PS or CO - exists, and what its effects are 
with regards to experimental probes. Four representative 
points along those two paths (see Fig^) were chosen and 
the corresponding spectral functions calculated. 

Figure|21a) shows A(k,ui) for the purely SC case (J = 
0) for /i = —1, leading to a uniform density (n) w 0.7. 
Fig. [21b) is for the case when the system presents lo- 
cal AF/SC coexistence (namely, both o.p.s simultane- 
ously nonzero at the same site) and Fig. |2Jc) for the 
pure AF phase. The red color indicates large spectral 
weight, whereas the blue one indicates very low inten- 
sity. In Fig. EJc), the AF gap can be clearly identified, 
together with the typical dispersion of the AF (upper 
branch), £k=±\/ e k + J 2 , which makes £k gapped ev- 
erywhere. This is in stark contrast to Fig. |2Ja), where 



the Fermi energy (Ep) close to (tt/2,tt/2), allowed by 
the symmetry of the pairing state. The "intermediate" 
state with local AF/SC coexistence is not drastically dif- 
ferent from the one with AF correlations only, and its 
resulting energy dispersion can be simply described by 
£k=±\/ (ek — m) 2 + J 2 + A 2 , once the parameter Ak is 
known. This conclusion is not supposed to change using 
the SO(3)-symmetric spin model. 

Similarly, along path 2 of Fig. ^ a point in the phase 
diagram with striped order was chosen, and the corre- 
sponding spectral density is given in Fig. Old). This 
result compares very well with previous calculations, 
(Ref. |H, Fig. 7): for instance, the system presents 
a FS crossing near (tt,0). Whereas the results from 
Fig-I3a),(c) refer to generally well-understood phases of 
the cuprate phase diagram, Figs.[21b),(d) are of relevance 
for the discussion related to the intermediate state, since 
they are both candidates for the intriguing phase in be- 
tween. 




FIG. 3: Experimental ARPES spectra for LSCO with x = 
and x — 0.03. Note the development of a (flat) second high- 
intensity branch near (n,0) and the emergence of a strongly 
dispersive signal at the Fermi level as the system is doped 
away from the half- filled insulator (reproduced from Ref. urn . 
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(0,0) (71,0) (71,71) (0,0) 

FIG. 4: (a) Dispersion £ k = ±^/(e k — y) 2 + A 2 . + J 2 for 
a perfect superconductor (SC) with J = 0, A = 0.3, a per- 
fect coexisting superconductor and antiferromagnet (SC+AF) 
with J = 1, A = 0.3, and a perfect AF with J = 1, A = 
as indicated. A k = — 2A(cos(fc :r ) — cos(fc y )) and e k is the 
usual dispersion of the free system, (b) Dispersion of a per- 
fect striped configuration. Lattice sizes are 12x12 and 8x8, 
respectively. 

For comparison. ARPES data from Ref. O 

for 

La2-xSr x Cu04 are reproduced in Fig. |3J For very low 
doping x = 0.03 (just inside the spin-glass insulating 
(SGI) phase) a flat band is observed close to -0.2eV in 
addition to a lower branch (energy ~-0.55eV), which is 
already present in the x = limit and therefore can 
be safely identified with the lower Hubbard band. As 
x is increased even further, the lower branch retains its 
energy position, but gradually loses its intensity until 
it is almost completely invisible after the onset of the 
SC phase at a;=0.06i&. In contrast, the second branch 
gains in intensity with doping, and also moves continu- 
ously closer toward the Fermi level; at the same time it 
starts to develop a coherence peak, which is clearly vis- 
ible for optimal doping. The main experimental result 
here, namely the existence of two branches near (tt,0), 
cannot be reproduced using spatially homogeneous models 
as demonstrated above. The cases of AF, SC and coex- 
isting AF+SC states all show only one branch below Ep 
nearby (ir, 0). This was already seen in Fig. |2Ja)-(c) for 
the MC data and is shown again in Fig. 01 a) for all those 
configurations (in all these cases the exact dispersion is 
known) . 

Tf _j_ • 12 :j i • 



(MC data) and Fig. 0{b) (perfect configuration of 
stripes), there will appear two branches near Ep, but 
the form of the dispersion is clearly different from the 
experimental data in Fig. [21 The investigation of A(k, w) 
for a spin-fcrmion model, related to Hup, but retain- 
ing the SO (3) spin symmetry, has been done carefully 
in Ref. |3^ Again, stripe phases were found for certain 
parameters and while in some cases the existence of two 
branches near (tt, 0) was reported, certainly there are no 
indications of "nodal" quasiparticles at (ir/2, ir/2). Then, 
stripes alone are not an answer to interpret the results 
of Yoshida et al. As a consequence, we conclude that 
neither local AF+SC coexistence nor stripes can fully 
account for the ARPES results in the low-doping limit 
and alternative explanations should be considered. 



Beyond the results already described, ARPES also pro- 
vides surprising insights/results for momenta other than 
(7r,0) (FigEJ. Along the Brillouin zone diagonal, a disper- 
sive band crossing Ep is found already in the SGI phase. 
The FS-like feature consists of a small arc centered at ~ 
(it/2, it/2); surprisingly, as more holes are added, this arc 
does not expand, but simply gains spectral weight. This 
increase in spectral intensity is roughly proportional to 
the amount of hole-doping for x<Q.l, although it grows 
more strongly thereafter. This observed increase in spec- 
tral weight is in relatively good agreement with the hole 
concentration rifj derived from Hall measurements and 
was interpreted as a confirmation of the hole transport 
picture. Below, however, we will provide a different ex- 
planation for this behavior. 



The aforementioned large gap (A«0.2cV) at (it, 0), 
together with the existence of the apparent gapless ex- 
citations around (tt/2, tt/2) is the essence of the PG 
problem. The shrinking of this gap and the concomi- 
tant appearance of a coherence peak has, for example, 
been interpreted as the evolution of a strongly coupled 
SC (at low doping) into a conventional BCS-SC at opti- 
mal doping. In this scenario, the large gap size directly 
reflects a large pairing scale, whereas the smallness of T c 
is attributed to the preponderance of phase-fluctuations 
in such a regime, which would outrule the existence of 
a phase-coherent SC condensate at higher temperatures. 
Alternatively, this gap may be regarded as the signal of 
a hidden order, which is not otherwise manifested. In 
other words, the relatively large excitation gap Apg is 
explained in terms of (i) a large SC gap Apq=Asc itself, 
or (ii) ApG=Asc+Aho, with a large, x-dependent hid- 
den order gap Aho whereas (iii) a mixed-state scenario, 
strongly influenced by disorder, leaves open the possibil- 
ity that it is the (local) chemical potential that deter- 
mines the PG physics. The precise role of fi in mixed- 
state phases needs to be examined further, but will not 
i_ _ .jj i i 




FIG. 5: Schematic representation of Sr doping. A chemical 
dopant (Sr) will not only disorder the nearest sites (blue color) 
in the CuCh-plane, but also neighboring ones, motivating the 
introduction of "plaquette"-like disorder configurations. 



B. Quenched Disorder I: Frozen Configurations 

Since calculations for A(k, u>) in the clean limit do not 
agree with ARPES measurements, we turn our atten- 
tion to a system with quenched disorder. The impact 
of quenched disorder is realized by tuning the coupling 
constants Ji and V\ in Eq.(0}. The spatial variations of 
these couplings is chosen in the following way (see Fig|5J): 
charge-depleted "plaquettes" that favor superconductiv- 
ity are placed on an AF background. The same procedure 
was followed in Ref.0, where more details can be found. 
The main point is that impurities will not only influence 
a single site, but, possibly due to poor screening in the 
proximity of the insulating phase, will change local po- 
tentials over a rather large area. The phase diagram of 
the clean model along path 1 and the corresponding dis- 
ordered case are reproduced in Fig.Elfor the benefit of the 
reader. Disorder has opened a region between the SC and 
AF phases where none of the competing order dominates 
and both regimes coexist in a spatially separated, mixed- 
phase state. This "glassy" state was discussed in detail 
in Ref . 0, where it was suggested that it leads to "colos- 
sal effects", in particular a giant proximity effect (GPE), 
which was recently observed in layered LSCO films 13 . 
The pronounced susceptibility of such mixed-phase states 
towards applied "small" perturbations is well-known and 
is, e.g., often regarded as the driving force behind "colos- 
sal magneto-resistance" in manganitesi^. 

To simplify the study and be able to access larger sys- 
tems, we will first consider a single SC cluster embed- 
ded in an AF background and also consider a fixed or 
"frozen" configuration of the classical fields (both AF 
and SC). Later, we will lift this restriction and perform 
a MC study. When a 12x12 SC region is placed on an 
AF background (total lattice size is 22x22), the resulting 
distribution of A(k, oj) is as shown in Fig.[7Ja). The con- 
tribution from the AF background is clearly distinguish- 
able from that of the SC island, since it is present even 
when the SC region is removed. The SC cluster induces 
a second "flat band" - quite typical for gapped systems 
- near E-p, along the (0,0) — > (tt, 0) direction. That this 
flat band is indeed produced by the SC island is verified 

1 J • - • - r _ • 1 1 j_ - o . . o /T^* IWI \ f-r. . i-T 
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FIG. 6: Results reproduced from Refll2lfor the benefit of the 
reader, (a) Phase diagram of model Eq. (|1J along Path 1 of 
Fig-Q showing the AF, SC and local coexistence regions (for 
spatially constant couplings), (b) Same as (a) when quenched 
disorder is added to the system in the form of SC plaquettes 
or impurities. In this case a region without long-range order 
appears. The temperature for PG formation is also indicated. 
Lattice size is 8x8 in both cases. 

(Fig. Ct) and finally for 5x5 (Fig. |7Ji), upon which this 
signal gradually decreases (the cases 9x9 and 11x11 give 
very similar results to 12x12 and are not shown.). The 
spectral intensity related to the surrounding AF "bath" 
concurrently decreases, in agreement with experimental 
observations^. 

The relative increase of the SC phase intensity also 
goes hand in hand with an increase in spectral intensity 
along the nodal direction, and the buildup of a FS around 
(7r/2, 7t/2). Figure |H1 shows a cut of A(k,uj) near Ep for 
the case depicted in Fig.EJa). There is considerable spec- 
tral weight near the FS for momenta close to (7r/2, 7t/2) 
only, suggesting that other parts of the FS arc gapped. 

Therefore, even the simplest possible mixed-phase 
state can qualitatively account for the observed ARPES 
data. It is also interesting to note that SC signals com- 
parably in strength with the ones stemming from the AF 
band, are only found for rather large SC blocks, encom- 
passing at least 20% space of the whole system. From 
this point of view, even in the strongly underdoped limit 
at £=0.03; the relative amount of the SC phase has to be 
quite substantial already. 

C. Quenched Disorder II: Monte Carlo Results 

The same system considered in the previous subsec- 
tion was evolved by the MC procedure explained in the 
introduction. However, the calculations can only be per- 
formed on smaller lattices. Results on a 10x10 lattice 
were obtained using a 4x4 region with couplings that 
favor superconductivity on a background that favors an- 
tiferromagnetism. As observed in the previous subsec- 
tion, the SC island produces the features S6en near , 
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(0,0) (tc,0) (%,%) (0,0) OcO) (7C.7C) (0,0) 

FIG. 7: Distribution of A(k,u>) for a single configuration 
of classical fields, corresponding to a SC region of size (a) 
12x12, (b) 8x8, (c) 7x7 or (d) 5x5 on a 22x22 lattice (i.e., 
30%, 15%, 10% or 5% SC respectively). Shown is E vs. k 
along (0,0) -» (tt,0) -» (tt.tt) -> (0,0). 




FIG. 8: Energy cut of A(k, w) close to Ef for a single con- 
figuration of classical fields, corresponding to a SC region of 
size 12x12 on a 22x22 lattice. The brightness of the col- 
ors indicate the intensity in the k x -k y plane for an energy 
cu = -0.241 



land, there is not enough resolution to be able to see 
the "flat band" described previously. This makes it un- 
likely for such structures to be seen in "exact" solutions 
of Hubbard-type Hamiltonians, which, due to a variety of 
reasons, are restricted to lattices as large or smaller than 
the ones considered here. As a further consequence, it is 
inevitable to conclude that the SC islands in the real ma- 
terial must be substantially large, of sizes approximately 
10 x 10 in lattice units which translates into 40Ax40A 

-1 i • 1 • _l 



D. Quenched Disorder III: Mean Field Theory 

The MC results presented above are also supported by 
more traditional approaches, namely the solution of the 
BdG equations. This self-consistent method allows for 
larger systems to be studied and is therefore much bet- 
ter suited to resolve the mostly subtle signals that are 
to be expected for the current investigation. In addi- 
tion, the tracking of the small FS as the doping level is 
changed certainly requires a lattice grid that allows for 
a reasonable amount of resolution in k-space. This, of 
course, cannot, at the current time, be achieved with the 
MC technique presented above, yet this does not mean 
this approach is without merit for the problem at hand. 
One has to keep in mind that the MC approach can be 
justified on the grounds that it is unbiased towards any 
particular state - at least, if performed in a correct fash- 
ion - whereas the self-consistent solution may fail in this 
regard in the "tricky" regime of small, but finite, dop- 
ing. The emergence of the stripe state in this case is well 
documented, but will only be realized for a sufficiently 
correct initial choice of the eigenfunctions. 

With these caveats in mind, it may be the best strat- 
egy to employ MC to establish the correct (clean) phase 
diagram, and then using this information to describe the 
emerging phases with the help of BdG. In this spirit, we 
present some results for the inhomogencous Hamiltonian 
i^HF below. Calculations were performed on systems up 
to A=32x32. As before, plaquettes were inserted to cre- 
ate areas of predominant AF and SC order, respectively. 
These plaquettes were randomly distributed throughout 
the lattice, and were allowed to overlap. Owing to the 
random distribution of plaquettes, clusters of low elec- 
tronic density are generated, their average size mainly 
determined by the the size of the underlying plaquettes 
(for the case of non-correlated plaquettes). We have cho- 
sen U\=h (AF sites) and Vj=-1 (SC sites), unless oth- 
erwise mentioned. With those values, it is possible to 
clearly separate the AF and SC signals, an important as- 
pect when studying not overly large systems, and when 
one is looking for subtle signals. Our main conclusions, 
however, are not supposed to change for other parameter 
values. 

The number of plaquettes grows linearily with the hole 
density, which ranges from £=0.03 to x~0.2 to mimick 
ARPES investigations. The total percentage of the area 
occupied by SC bonds is roughly proportional to the hole 
density; however, given the phase diagram in Figd there 
is a strong connection between the area, asc, occupied by 
the SC regions, and the area occupied by the AF phase, 
aAF=l-asc: 

1 - x = 1 * <2af + n sc * a sc , (10) 

where nsc is the SC density (the AF density is nAF~l), 
which leads to the following expression 

ten = . (11) 



Eq. can be interpreted in two ways: for a given 

value of x, cither a desired value asc defines rise- or vice 
versa. Here, our parameters (i.e. n\) are tuned so that 
ngc^O.75, the value upon where the system becomes a 
homogeneous SC according to FigQ] 

In Fig|5]we show typical configurations for which spec- 
tral functions were calculated, each corresponding to dif- 
ferent densities, (n)=0.97, 0.94, 0.87 and 0.79, respec- 
tively. The uppermost row ((a)-(a3)) displays the envi- 
ronment created by the plaquettes, with the red color 
favoring AF and the blue one the SC phase. In the sec- 
ond row ((b)-(b3)), the corresponding SC gap amplitudes 
are shown; yellow denotes large Aj's (which are approx- 
imately in line with the ones found for a pure system at 
density nsc) and ever darker colors mean an ever smaller 
value of Aj . The bottom row plots the corresponding AF 
o.p., which, given our choice of plaquettes, results in what 
is essentially a mirror image of the second row. Here, blue 
denotes strong local AF amplitude mi, whereas lighter 
colors stand for very weak or absent magnetic ordering 
((c)-(c3)). What the second and third row in FigEjl show- 
is that there are essentially three different regions in such 
a mixed-state: one that is a pure SC with no discernible 
AF amplitude, one where both amplitudes assume finite 
values (those may vary quite distinctively) and a third 
one that is solely AF. Those intermediate regions appear 
even though the underlying plaquettes are either purely 
AF or SC (which probably is an oversimplification to be- 
gin with), and owe their existence to boundary effects. 

The spectral functions related to the mixed AF/SC 
states such as in FigEDare shown in Fig^| in addition 
with those for the pure AF (x=0) case (Fig llOf a)) and the 
pure SC (nsc=0.76, Fig lf Of f)): the former has the usual 
dispersion, characterized by a flat band in the neighbor- 
hood of the X point, plus the shadow bands near (ir, 
7r). In Fig |10f a1-(Tl. the yellow color stands for large 
spectral weight (with respect to the maximum of A(k, uj) 
for each density), with ever darker colors describing ever 
smaller intensities. With the introduction of SC pla- 
quettes to the AF background fFig lTUT al— >(b)). spectral 
intensity is accumulating close to Ep at momenta (n, 0) 
and at E-p for ~(7r/2, 7r/2), both of which were gapped 
at x=0. At the same time, the dispersion arising from 
the lower Hubbard band remains clearly identifiable and 
almost unchanged in position at energy u>AF^-2t. There- 
fore, in the neighborhood of the X point two distinctive 
branches are emerging for phase- separated systems such 
as in FigO and both the SC and the AF branch can 
be clearly identified by comparison with the respective 
"clean" regime (x=Q, x~0.2). Their relative intensities 
depend on the ratio ciaf/osC: with the AF branch be- 
coming increasingly faint as this ratio tends towards 0. 
The observed features are very broad for weak doping 
for any wavevector close to E-p fFig llOf bhfc')'). charac- 
teristic of poorly defined quasiparticles with very small 
residue Z . This residue, however, appears to be small be- 
cause of the disordered nature of the ground state, rather 




FIG. 9: Typical quenched disorder configurations on a 32x32 
lattice at different electronic densities: the top row (a)-(a3) 
displays the local environments (red=AF, blue=SC), whereas 
(b)-(b3) shows the emerging local SC amplitudes, with the 
yellow color standing for the highest values (black for Ai=0). 
The bottom row shows the AF amplitudes. Blue equals strong 
AF ordering, whereas ever lighter colors identify a weak AF 
order parameter. 
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FIG. 10: AQi, ui) vs. k calculated for a series of densities 
on a 32x32 cluster, starting from x=0.0 (100% AF) in (a), 
to x=0.03 (22% SC) (b), x=0.06 (43% SC) (c), x=0.13 (70% 
SC)(d), x=0.21 (84% SC) (e), and the 1=0.24 (100% SC) 
(f). The two branches visible near (ir,0), belonging to the 
SC and AF region, respectively, can be clearly identified for 
the intermediate states, which suggests a simple mixed-state 
interpretation for ARPES experiments. 
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FIG. 12: Integrated spectral weight Jg, for two different 
integration windows (in k-space, with frequency cutoff loo—- 
0.3) as a function of the hole density x. The downturn of 7^ 
at large doping is a finite-size effect. 



FIG. 11: The spectral intensities, integrated over a small 
shell below the Fermi level as described in the text, at dif- 
ferent densities (n)=0.97 (a), (n)=0.90 (b), (n}=0.87 (c) and 
(n)=0.79 (d), and using a 32x32 lattice. The FS arc devel- 
oping around ~ (0.457T, 0.457r) can be readily identified. It 
increases in intensity only as holes (plaquettes) are added to 
the system, whereas the arc's length remains essentially con- 
stant. This is to be compared with ARPES measurements in 
underdoped LSCO±i. 



charge carriers - quite an important distinction! For ex- 
ample, as one moves closer towards the homogeneous 
state f Fig HoT o) — » ( d) ) . the intensity distribution at (tt, 
0) becomes much sharper and reminiscent of a quasipar- 
ticle peak, without changing the interactions at all. All 
those observations are in agreement with Fig0 (although 
those data are much better defined, presumably because 
the o.p.'s were assumed constant for each region), as well 
as with experimental evidence, which has uncovered a 
very similar behavior in a series of underdoped to opti- 
mally doped LSCOia. 

Similar agreement between theory and experiment is 
found around the (n/2, ir/2) point, where significant 
spectral weight akin to a quasiparticle peak is found right 
at the Fermi level, and this weight continously increases 
as x grows. This can be seen in Fig llOl but is more 
obvious in Fig^] which shows the spectral intensities, 
integrated between the chemical potential and a cutoff 
WD=-0.3t, at different densities. Note that these data 
are calibrated with respect to the maximum intensity as 
found in Fig lllf d). The emergence of a FS 'arc', cen- 
tered in the neighborhood of T, is obvious. With increas- 
ing doping it is the arc intensity that increases, whereas 
the arc size barely changes (compare, e.g., Fig lTlT c).(d')'). 
Both observations mirror ARPES results and suggest, 
in conjunction with these data here, a simple picture of 
the underdoped cuprates involving AF and SC cluster- 



calculated the integrated spectral intensity in the neigh- 
borhood of the FS at different doping levels, Jg: 

Jg = f A(k,u)d(jdk, (12) 

J-\u>n\ JO 

where □ denotes a window ("window 1(2)" is a 5 2 (7 2 ) 
square) centered around the FS arc midpoint. Using data 
gathered from ARPES experiments (however at wd~0), 
it was shown that this quantity increases approximately 
lincarily with doping 11 , which is also what arises from the 
mean- field calculation (Fig ll2fl . This increase is almost 
entirely associated with increased spectral weight around 
the nodal direction, as is also clear from Fig(nj Then, 
the most natural explanation for the observed behavior 
of A(k, Lu^E-p; x) is that it simply reflects the increasing 
amount of age rather than the change in particle density 
(although both arc related). A linear relationship be- 
tween Ijj and the particle density as such is not obvious 
(or even fulfilled) for a clean SC system, but it is natu- 
rally explained in a mixed-state picture. For the perfect 
SC (nsc~0-76), the FS is even more clearly defined as in 
Fig llir d'): in this case, the central peaks remain and gain 
in intensity, whereas the small incoherent intensity in the 
rest of the fc^-A^-plane vanishes. Again, this "sharpen- 
ing" of quasiparticle peaks as the system becomes more 
homogeneous - and quasiparticles better defined - has 
been observed in ARPESit 

We also want to mention here that the experimen- 
tal observation - a doping-independent size of the FS 
arc - also contradicts other popular proposals, which as- 
sociate the underdoped regime with one of strong pair 
attraction and significant phase fluctuations: In those 
scenarios, V increases as x— >0, but this would entail a 
FS arc that at the same time is shrinking in size. The 
simple fact that the size stays constant suggest a rela- 
tively doping-independent value of A and, thus, V. The 
same is true for scenarios, which rely on additional, non- 
zero order-parameters (gaps) dominating the underdoped 
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this additional order becomes stronger at lower doping, 
in contrast to what is seen experimentally. 

To lend further insight to our analysis, we have also at- 
tempted to a series of different calculations at fixed den- 
sity (n)=0.87, but somewhat deviating from the above 
scenarios: (i) first, a model where the random SC config- 
uration as shown in Fig[!Ja2) is replaced by one with a 
single cluster, occupying the same overall area, in order 
to better understand the effects of randomness, cluster- 
size and -shape; (ii) second, a model of charge-depleted 
plaquettes with V=0, i.e. a situation of hole-rich and 
hole-poor phases coexisting, but without the former one 
possessing SC order, and finally (iii) a model, where pla- 
quettes are distributed in an ordered fashion, forming a 
supcrlatticc of 16 plaquettes, each covering a block of 
6x6 sites on the usual 32x32 lattice. A(k, u)'s result- 
ing from those configurations are shown in Fig ll3f bWdh 
In the case of (i), the replacement of the random con- 
figuration by a "superblock" leads to a narrowing of the 
spectral functions, with more clearly defined peaks, but 
otherwise leaves the overall features such as the existence 
of the two-branches and the FS positions, unperturbed. 
One consequence, however, is the narrowing of the gap 
distribution function P(A) as the randomness is reduced, 
reflecting the smaller number of possible environments. 
For the present case, P(A) consists of a central peak from 
data well inside the cluster and a smaller peak from the 
SC/AF boundaries, whereas for the totally random case 
P(A) can be approximately described by a Gaussian dis- 
tribution, with a multitude of A's. This might be of 
importance, since exactly such a Gaussian distribution 
has been reported in STM experiments in underdoped 
BSCCO. Then, randomly located SC regions are the best 
picture to describe the experiments. On the other hand, 
quite drastic changes occur if the plaquettes are chosen 
as charge-depleted-only regions, as shown in Fig llbf c): 
in such a case a single band with a clear FS crossing 
for k~(7r,0) emerges, and such a feature is certainly not 
observed in experiments. Therefore, the inclusion of a 
pairing term even in the low-doping limit seems neces- 
sary to correctly reproduce the experimental features. 

The case of the ordered plaquettes fFig lCT dl) was mo- 
tivated by claims that the STM results in underdoped 
(but SC) BSCCO should be best interpreted in terms 
of charge-ordering^ the superlattice structure chosen is 
the simplest realization of such a CO state, and may 
help elucidate the ARPES measurements of such a com- 
plicated, yet highly ordered, phase. Although the overall 
features are largely unchanged when compared to the 
random cluster situation, the broad peaks observed for 
the random structure split up into several subpeaks, i.e. 
for some momenta there may be even more than just two 
solutions. This is even clearer if more traditional repre- 
sentations of A(k, lo) are chosen, not shown here. Al- 
though those multiple-peak features may be too weak to 
be observed in current ARPES experiments, and there- 
fore such a state cannot a priori be ruled out, there is no 
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FIG. 13: yl(k,aj) for doping rate £=0.13: (a) same as 
Fig llOf c'l. (b) has one single SC-favoring cluster in the middle, 
(c) Vi=0 everywhere, even in charge-depleted regions, and (d) 
with SC plaquettes forming a super-lattice. The SC gap dis- 
tributions -P(A) are shown in (e) for the random model (a) 
(thick black line), the single cluster (b) (dashed line) and the 
superlattice structure in (c) (red line). In the first case, the 
gaps have a Gaussian-like distribution around the maximum, 
with a broad tail, whereas in case (b) there are essentially two 
peaks only, (c) has a multitude of clearly defined peaks. 



least according to ARPES investigations. 



IV. DOS AND T* 

Calculations of the total DOS for the model with 
quenched disorder is shown in Fig. 114b . This is at a 
point in the phase diagram, Fig. HJb) with 6 impurities, 
i.e. without long-range order either in the AF or the SC 
sector. The DOS clearly shows a PG that disappears 
at a temperature scale T*. The qualitative physics de- 
scribed here has already been investigated by the authors 
in Ref. an d will not be repeated. 

It is worth noting that in this model the PG is caused 
by short-range order of either SC or magnetic variables. 
In other words, if the temperature for short order forma- 
tion is denoted by T* c and T* { for e?-wave superconductiv- 
ity and antiferromagnetism, respectively, then the tem- 
perature for PG formation is roughly T*^max(T s * c , T* f ). 
This explain the dependence of T* with doping as ob- 
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short- (and possibly long-range) range AF order and no 
SC order, therefore, T* ~ T* { for low doping. As doping 
increases, the system becomes less and less AF which im- 
plies a decrease in T* f and consequently in T*. However, 
as the carrier concentration is increased further, the sys- 
tems starts to present SC order, at short distance first 
and then at long distances. Then, T* stops decreasing 
and stays constant or even increases with carrier concen- 
tration near the optimal doping region. 

The BdG equations are less influenced by finite-size ef- 
fects and, thus, allow for a better understanding of the 
DOS of a mixed AF/SC-statc. Although we have only 
performed TwO-calculations only, some interesting infor- 
mation can be extracted from these solutions as discussed 
below. As the system is doped away from half-filling, the 
SC regions appear as mid-gap states between the original 
Hubbard bands (Fig^Jb)), which are gradually filled as 
doping is increased, at the expense of the original states. 
The associated peaks in N(u>) gain in strength, whereas 
the energy states belonging to the AF sites gradually 
fade away, giving way to the usual DOS of a dSC. From 
the most naive point of view, this two-peak structure in 
the DOS can be associated with the PG, wherein it is 
the disappearance of the AF peak (at higher energies) - 
which should be happening at ever lower T's as x is in- 
creased - that determines T*(x). But there are also some 
other subtle issues that need to be considered following 
Fig ll-lf bl: the SC peak seems to be travelling towards 
higher energy, rather than towards E-p. This can already 
be seen in the corresponding ARPES data, but is much 
more obvious in the DOS. This is in disagreement with 
data from ARPES, which clearly show the SC band mov- 
ing closer to Ep as x is increased. In this framework, 
the PG is dictated by the energy position of the charge- 
depleted phase in the presence of the half-filled insulator. 
This is an important, yet subtle issue since the energies 
involved are rather small; yet it should be explored in 
further considerations of the mixed-state picture. Sum- 
marizing, although the description of PG physics is not 
perfect, at least the key features (such as its existence) 
are neatly reproduced by our mixed-phase model. 



V. LOCAL DOS 

Consider again a system with inhomogeneous cou- 
plings Ji and V\ in a such a way as to produce AF 
and SC regions in the sample. We calculate Ni(w) for 
all sites, displaying them in "linear order"— in Fig. I15t 
1161 for various plaquettes concentrations as indicated. 
The results show that there are two types of curves or 
"modes", and that each site contributes to a different 
"mode" of the LDOS. For example, all AF sites present 
a clear gap (Fig. H5k ). whereas SC sites have the equiv- 
alent of a d— wave gap (on a finite-size system). This 
seems in agreement with STM measurements^ in un- 
derdoped Bi 2 Sr2CaCu208+<5- For example, Figure 3 of 



25 



20- 



3 



15> 



10- 



(a) 



/? //• 

ft'//: 

/A. 



-T= 


=0.02 


-T= 


=0.04 


- T= 


=0.06 


-T= 


=0.08 


-T= 


=0.10 


T= 


=0.14 


-T= 


=0.24 



J L 



-0.8 -0.4 





CO 



0.4 0.8 



— <n>=l 

— <n>=0.97 
-- <n>=0.94 

— <n>=0.87 
<n>=0.76 




FIG. 14: (a) DOS near the Fermi energy (0) for the situation 
of Fig. |S| where there are 6 plaquettes on an 8x8 lattice. The 
disappearence of the PG can be seen as the temperature is 
increased, (b) DOS calculated at different densities by mean- 
field theory at T=0 for a 32x32 lattice for the disordered 
model. The clean case ((n)=0.76) was calculated for a 40x40 
lattice to minimize finite-size effects. 



on the sample vs. the bias, indicating two types of regions 
- similarly to what is observed in Fig. 115b and Fig. 116b 
of the present work, which represent intermediate states 
which do not present a dominant global order parameter 
(as is the case for Figs ll5f a'). ll6f a'). Certainly these data 
still suffer from finite-size problems, but the basic issues 
are captured nevertheless. Also, note that the transition 
from the SC to the AF regions is gradual with an inter- 
mediate phase that has a small SC-like gap, but lacks the 
pronounced (coherence) peak of the SC state. These re- 
sults favor the idea that "exotic" phases are not needed 
to explain the nature of undcrdoped cupratcs, but in- 
stead that quenched disorder creates regions of local SC 
or AF order when those two phases strongly compete. 

The distribution of intensities of the LDOS for each 
frequency range is shown in Fig. 1171 The AF, SC and 
metallic (non-SC) contributions are present with differ- 
ent intensities depending on the concentration of SC pla- 
quettes as shown. A metallic, but not SC, phase appears 
due to the value of the chemical potential used around 
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FIG. 15: Local DOS Ni(uj) vs. uj — fi on an 8x8 lattice with 
1 and 6 "plaquettes" , respectively, plotted in "linear order". 
The color convention is as follows: red indicates AF sites; 
orange, SC sites and cyan, metallic sites. The inset shows the 
ei-wave gap distribution, |Ai, x |(cos((^i) — cos(<^>i +x )) in each 
case. 



FIG. 16: Ni(u)) on an 8x8 lattice with 20 and 24 "plaque- 
ttes" , respectively plotted in "linear order" . Color convention 
as in Fig |15l Just as in Fig |15l the datasets for each site i are 
offset by a small amount from each other to allow for better 
visibility. 



data presented here, but instead provides a more realistic 
separation between SC and AF regions. With a simpler 
model for the plaquettes, the metallic peak would not be 
present. 

This is even more obvious in the mean-field approxima- 
tion, where sites with local AF or SC order can be easily 
identified. Travelling along certain cuts in a mixed-phase 
system, which show either AF or SC order, and measur- 
ing Nj(ui) along this path (see Fig H7T b).(c)). it is re- 
vealed that the signals are vastly different in either case. 
Whereas the AF sites (derived from the value of mi have 
a clear and well-defined gap, the SC sites resemble those 
of a d-wave SC. From these results it is also clear that 
there is not just one SC gap (which may, for example, 
be extracted from the position of the first main peak), 
but rather a distribution of gap values, presumably re- 
flecting the random nature of the samples, as depicted 
in Fig|5Jb)-(b3). As the hole doping level is increased, 
the relative amount of SC data taken "in the bulk" in- 



Those data, for example, should be compared with those 
for Cai.92Nao.o8Cu02Cl2^, where similar features were 
observed; sites with a large excitation gap were found 
to be insulating, whereas the other ones were metallic. 
As the number of charge-carriers is modified, the relative 
amount of those two phases changes, whereas the dl/dV- 
signals themselves remain fairly independent of the dop- 
ing level. 



VI. CONCLUSIONS 

In a recent publication, a phcnomcnological model for 
cuprates was introduced^. This model includes itiner- 
ant fermions coupled to classical degrees of freedom that 
represent the AF and SC order parameters. The model 
can be studied integrating out the fermions and Monte 
Carlo simulating the classical fields. The simple charac- 
teristics of this model allow to investigate the crossover 
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FIG. 17: Distribution of the local DOS vs. w - ^ on an 
8x8 lattice with the concentration of "plaquettes" shown and 
parameters corresponding to Figs. 1151161 NI stands for the 
number of impurities, (b) shows Ni(u) for sites, which are 
predominantly AF, whereas (c) those for SC ordered sites. 
Those data were taken on a disordered lattice with 28x28 
sites, and about 60% SC. The parameters were f/i=4 for the 
AF regions, and Vi=-1 for the SC sites. 



near the Fermi energy for undcrdoped compounds are 
explained within the context of our model as induced by 
the SC regions. 

The calculation of the DOS revealed a PG for the 
regions of the phase diagram where no long range or- 
der dominates, although local order exists. The disap- 
pearance of the PG with temperature gave an estima- 
tion of the temperature scale, T*. Moreover, the LDOS 
calculated in the inhomogeneous system presents two 
types of curves or "modes" , corresponding to AF and 
SC clusters. This observation is in agreement with STM 
measurements^. 

In summary, a simple model is able to capture several 
features found experimentally in high temperature super- 
conductors. In this theory, the glassy state is forme d by 
patches of both phases, as already discussed in ReflTJ. 
and it can present giant responses to small external per- 
turbations. In this work, it was shown that ARPES, DOS 
and LDOS experimental information in the undcrdoped 
regime can also be rationalized using the same simple the- 
ory. Our model does not address directly the important 
issue of the origin of pairing in the SC state, but can de- 
scribe phenomcnologically the AF vs. SC competition. 
The results support the view that undcrdoped cupratc 
superconductors are inhomogeneous at the nanoscalc, 
and that only well-known competing states (AF and SC, 
perhaps complemented by stripes) are needed to under- 
stand this very mysterious regime of the cuprates. 



porated, and regardless of whether homogeneous or inho- 
mogeneous states emerge as ground states. The simpli- 
fied character of the model, as compared with the much 
more difficult to study Hubbard model, allows for numer- 
ical studies at any electronic density and temperature, 
and the evaluation of dynamical properties as well. Due 
to present day limitations in the analysis of many-body 
problems, the complex physics of transition metal oxides 
at nanometer length scales can only be captured with phe- 
nomenological models, such as those recently discussed by 
Alvarez et alH. 

In the present publication, the previous effort has been 
extended to the analysis of photoemission one-particle 
spectral functions and the local DOS, comparing theory 
with experiments. It has been observed that without 
quenched disorder (clean limit) the model presents var- 
ious phases: SC, AF, local coexistence of SC and AF 
orders, and striped states. However, the spectral density 
calculated for all these states does not reproduce the ex- 
perimental measurements for cuprate superconductors as 
reported in Ref. 0. To make progress in the theoretical 
description of experiments, quenched disorder was added 
to the system, inducing regions of SC and AF order, in- 
terpolating between the two fairly uniform AF and SC 
states. In this case, the ARPES spectral weight presents 
two branches, one induced by the AF background and 
the other by the SC regions or islands. This spectral 
weight shows clear similarities with the experimental ob- 

j_* i j_i r - a • _i i 
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APPENDIX: DIAGONALIZATION PROCEDURE 

To diagonalizc Eq.Q) a modified Bogoliubov 
transformation^ needs to be applied. After some 
algebra it can be shown that Eq. (0J becomes: 

N 1 

H = Y Y XX+^i^iW^'+^jHlt^'T 

n,n' — 1 a.a'—O (ij) 
N 1 

~ Y Y Y b n+o<N( i ) H bj a > b n'+a'N{jhn>ni l + 
n^n' — l a,a'— (ij) 

N i i 
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where W and H * arc the 2N x 2N matrices given by 

K (A) is a N x TV matrix, such that JVy = —t (Ay = 
Ai, Q ) only if i and j are n.n. sites and otherwise, 
iy = JiSfdy and M„ = E U ^(i)(^y - Jy)Mj) + 
a n+Ar(i)(^ij _ ) a n+N (j) ■ The c-numbcrs a„(i) (6„(i)) 
are chosen, so that they diagonalize {H+), i.e.: 

'E a>n+aN {i)Hj a ^ a i an' +a' N {}) = E n+aN 5 n + a N : „< +a i jy ■ 
ij 

(A.4) 

Then, the total energy can be written as 

n=N 

Etotal = Yl^ E lf^ E D +Ei J^ E n)) 
n=l 

N 1 1 

+ ^(M, l - J B.,i) + -^-|A iQ | 2 , (A.5) 

n=l i.a 1 

where f(x) = 1/(1 + c x ) is the Fermi function. The sum 
is only over the N largest eigenvalues of W, {i?T}i< n <j\r, 
and the N largest eigenvalues of H^, {£!^}i<„<jv- This 
expression involves the eigenvalues of both matrices. Al- 
ternatively, it is possible to express -Etotal in terms of the 
eigenvalues of only one matrix, for example, which is 
the more efficient way for the MC simulation. To under- 
stand that, first note that the eigenvalues of have the 
opposite signs as those of W. Let : 




(A.6) 



where / is the N xN identity matrix. Then S = S^ 1 and 

= —SH^-S -1 . Therefore, if \v) is an eigenvector of 
with eigenvalue E^, then S" 1 ^} is an eigenvector of 
with eigenvalue —E\. This proves that the eigenvalues of 
H * and H* arc the opposites of one another. From this 
discussion it also follows immediately that if a is an eigen- 
vector of H\ then b defined by b n = a n+ N, b n +N = a>n 
Vn G [1, iV] is an eigenvector of so the eigenvec- 
tors of one matrix can be obtained from the eigenvectors 
of the other. Now, let {EI}n+i<u<2N be the TV low- 
est eigenvalues of and similarly {£^}jv+i<n<2JV the 
N lowest eigenvalues of HK Therefore, E^ — — E^ +N 
Vn G [1, AT]. Then, the second term in Eq. I|A.5() becomes 

- E l+Nf(-@ E i+N) and usin S that f(~ x ) = 1 - f( x ) : 

n=2N N 

^totai = e E if(P E i) + E M » + 2 £ ^ |AiQ|2 ' 

n— 1 n— 1 i,a 



and this expression was used in the MC evolution of the 
system. 

Most observables are calculated by replacing the elec- 
tron operator c- la by Bogoliubov operators via Eq. JSJ. 
For example, the number of particles, N e , is given by the 
average of c[ a c\ a and in terms of a n (i) its expression 
is: 

n=2JV N 
N e = E l«n| 2 + 2E(l a ™| 2 -| a «+^| 2 )(/(^i) + /(^))' 

71 — 1 71 — 1 

(A.8) 

where we have used the abreviation |a„| 2 = |a„(i)| 2 . 
Note that unlike the standard Bogoliubov expression for 
the number of electrons, the second term in Eq. IA.8I) 
can contribute even at T = due to the fact that E^ 
can be negative for J finite. Moreover, unlike standard 
spin-fermion models, such as those for manganites and 
cupratesjiS the number of electrons depends not only on 
the eigenvalues of the one-particle sector, but also on the 
eigenvectors. 

As a particular case consider A(r,t), which is defined 
by the expression: 

A(r, t) =< E 4(*)<W(°) + H.c. >, (A.9) 
l 

where (...) denotes thermal averaging. Applying the 
modified BdG transformation, Eq. (JSJ, Eq. (|A.9|I is cal- 
culated using: 

A{v^) = Y J X n {T)8{u-El) + Y n {v)5{w + El), (A.10) 

n 

where 

X n (r) =E«:aK(l + r), (A.ll) 
1 

and a similar expression is valid for Y n . Eq. (|A.10() can 
be Fourier-transformed to obtain A(k, to), but it is faster 
to do that after performing the average, and that route 
has been followed in the present work. 

In STM experiments, one typically measures the 
change of a (local) tunneling current dl/dV, a quantity, 
which - under suitable assumptions - is proportional to 
N(i, uj) and thus allows for a direct mapping of the local 
electronic states. N(i,u>) is but the Fourier-transform of 
A(k, u>) and can be just as easily evaluated: 

N 

N(i t w) = E \a n (i)\ 2 S{u - E n ) + \a n+N (i)\ 2 5(u + E n ). 

71 = 1 

(A.12) 
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